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r— i ABSTRACT 

Taking advantage of strong gravitational lensing, we have recorded the rest-frame ultraviolet 
\ spectrum of the z = 2.38115 galaxy 'The Cosmic Horseshoe' (Jl 148+1930) at higher reso- 

• lution and signal-to-noise ratio than is currently feasible for unlensed galaxies at z = 2 — 3. 

^ ■ With a star-formation rate of SFR ~ 1OOM yr -1 , dynamical mass M vir ~ 1 x 1O 1O M , 

half-solar metallicity, and moderate reddening E(B — V) = 0.15, the Cosmic Horseshoe is a 
Q ■ good example of the population of galaxies responsible for most of the star-formation activity 

' at these redshifts. 

From the analysis of stellar spectral features we conclude that a continuous mode of star 
formation with a Salpeter slope for stars in the mass range 5-100Mq gives a good represen- 
tation of the UV spectrum, ruling out significant departures from a 'standard' IMF. Generally, 
— . we find good agreement between the values of metallicity deduced from stellar and nebu- 

1 lar tracers. Interstellar absorption is present over a velocity range Av ~ 1000km s _1 , from 

[ —800 to +250km s _1 relative to the stars and their Hll regions, but we still lack a model 

relating the kinematic structure of the gas to its location within the galaxy. There is evidence 
to suggest that the outflowing interstellar gas is patchy, covering only ~ 60% of the UV stellar 
continuum. 

\Q ■ The Lya line shares many of the characteristics of the so-called Lya emitters; its double- 

' peaked profile can be reproduced by models of Lya photons resonantly scattered by an ex- 

0^ . panding shell of gas and dust, with ~ 10-15% of the photons escaping the galaxy. Many 

of the physical properties of the Cosmic Horseshoe are similar to those of the only other 
| galaxy at z = 2-3 studied in comparable detail up to now: MS 1512-cB58. The fact that these 

• *h . two galaxies have drastically different Lya lines may be due simply to orientation effects, or 

/\ ' differences in the covering factor of outflowing gas, and cautions against classifying high-z 

\ galaxies only on the basis of spectral features, such as Lya, whose appearance can be affected 

by a variety of different parameters. 

Key words: cosmology: observations — galaxies: evolution — galaxies: starburst — galax- 
ies: individual (Cosmic Horseshoe) 



1 INTRODUCTION 

In the thirteen years that have elapsed since the term "Lyman break 
galaxy" (or LBG) was coined, samples of galaxies at z — 2~A 
have increased a thousand-fold. Large surveys, selecting galaxies 
not only via the break at 912 A (due mostly to intergalactic H I ab- 
sorption), but also via Lya emission, a variety of colour criteria, 
and sub-mm emission, have given us a broad view of the galaxy 
population during the epoch when star-formation activity was at its 
peak in the history of our Universe (see Pettini et al. 2008 for a 
brief review). 
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However, typical galaxies at z ~ 3 are faint, with fiducial 
7£-band magnitude = 24.4 (Steidel et al. 1999; Reddy et al. 
2008), so that, even with the most efficient optical and infrared (IR) 
telescopes and spectrographs currently available, it is only possi- 
ble to record their spectra with limited resolution and signal-to- 
noise ratio. Until the next generation of 30+ m optical/IR telescopes 
comes into operation, the only strategy at our disposal for studying 
in detail the internal properties of normal galaxies at these redshifts 
is to take advantage of fortuitous alignments with foreground mass 
concentrations (massive galaxies or galaxy clusters) which can pro- 
vide an order-of-magnitude boost of the flux reaching the Earth 
through gravitational lensing. 

For several years, the archetypal object for this kind of study 
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has been the Lyman break galaxy MS 1512-cB58 (or cB58 for 
short; Yee et al. 1996). The magnification by a factor of ~ 30 (Seitz 
et al. 1998) provided by an intervening galaxy cluster at z = 0.37 
brings this L* , z — 2.7276, galaxy within reach of high resolution 
spectroscopy from the ground and has made it the target of exten- 
sive observations at wavelengths from the visible to the mm range 
(Pettini et al. 2000, 2002; Savaglio, Panagia, & Padovani 2002; 
Teplitz et al. 2000; Siana et al. 2008; Baker et al. 2001, 2004; Saw- 
icki 2001). These studies provided unprecedented clear views of 
the interstellar medium, young stars, star-formation history, dust, 
metallicity, kinematics, and many other physical properties of a 
'typical' LBG at z ~ 3. 

But how typical is cB58? We know from the work of Shapley 
et al. (2003) that its ultraviolet (UV) continuum is more reddened, 
and its interstellar absorption lines are stronger, than in an 'aver- 
age' LBG. Furthermore, its young age of only 10-20 Myr (Elling- 
son et al. 1996; Siana et al. 2008) is at the lower end of the wide 
range of values deduced by fitting the spectral energy distributions 
of galaxies at z — 2-3 (Papovich, Dickinson, & Ferguson 2001; 
Shapley et al. 2001; Erb et al. 2006b), although it is likely that an 
older stellar population is present too. Clearly, it is important to 
study in detail several other similarly bright sources and establish 
the range of properties of galaxies at a given cosmic epoch, rather 
than naively assuming that any particular galaxy is representative 
of a whole population. Fortunately, the large area of sky surveyed 
by the Sloan Digital Sky Survey (SDSS) has led to the recent identi- 
fication of many other examples of highly magnified galaxies (see, 
for instance, Belokurov et al. 2009; Kubo et al. 2009, and refer- 
ences therein). These and other discoveries have spurred a num- 
ber of follow-up studies (e.g. Smail et al. 2007; Stark et al. 2008; 
Finkelstein et al. 2009; Hainline et al. 2009; Yuan & Kewley 2009). 

With these new samples of objects it will be possible to make 
progress on a variety of outstanding questions which cannot be sat- 
isfactorily addressed with lower resolution data. In particular: 

(i) To what extent do different metallicity indicators, based on 
stellar photospheric lines, wind lines from the most massive stars, 
interstellar absorption lines, and emission lines from H II regions, 
give consistent answers? This is an important issue, not only to ex- 
plore the degree of chemical homogeneity of galaxies undergoing 
rapid star-formation, but also to clarify the sources of the system- 
atic offsets between different H II region metallicity calibrators, and 
possibly bring them into better internal agreement (see, for exam- 
ple, the discussions by Pettini 2006, and Kewley & Ellison 2008). 

(ii) What do the profiles of the interstellar absorption lines tell 
us about large-scale outflows in actively star-forming galaxies and 
about the inflow of gas fueling star formation? Can we see any 
evidence of the leakage of hydrogen ionizing photons from the sites 
of star formation into the intergalactic medium (IGM) which has 
proved so difficult to detect directly (e.g. Shapley et al. 2006; Iwata 
et al. 2009) and yet seems to be required by a number of indirect 
lines of evidence (e.g. Faucher-Giguere et al. 2008)? 

(iii) Can we place limits on possible variations of the stellar 
initial mass function (IMF) at z — 2-3? Claims to this effect have 
certainly been put forward (e.g. Wilkins et al. 2008 and references 
therein), but the evidence is still controversial (e.g. Reddy & Steidel 
2009). 

(iv) What are the factors affecting the wide variety of spec- 
tral morphologies of the Lya line, from strong narrow emission 
to damped absorption (e.g. Mas-Hesse et al. 2003; Verhamme, 
Schaerer, & Maselli 2006), and what can we learn about the re- 
lationship of the so called 'Lya emitters' to the more general LBG 
population? 




Figure 1. /2-band image of the Cosmic Horseshoe obtained with FORS2 on 
the "Very Large Telescope of the European Southern Observatory (courtesy 
of L. J. King). North is up and East is to the left. Superposed on the image 
is the 1.0 X 20 arcsec entrance slit of ESI used for the observations reported 
here. 

In this paper we address these and other related questions with 
new observations of the rest-frame UV spectrum of a recently iden- 
tified highly magnified LBG, named 'The Cosmic Horseshoe' by 
its discoverers (Belokurov et al. 2007). The paper is organised as 
follows. In Section 2, we summarize the known properties of the 
Cosmic Horseshoe and in Section 3 we provide details of our ob- 
servations and data reduction. In Section 4, we consider the infor- 
mation provided by the stellar UV spectrum of this galaxy, con- 
cerning in particular its systemic redshift, massive star IMF, and 
metallicity of the OB stars. Section 5 deals with the resolved pro- 
files of the interstellar absorption lines, while Sections 6 and 7 fo- 
cus on the UV nebular emission lines covered by our spectrum, Lya 
and C III] A1909 respectively. In Section 8, we bring these different 
strands together in the light of the questions raised above. Finally, 
we summarize our main conclusions in Section 9. Throughout the 
paper, we adopt a cosmology with Q,m = 0.3, J1a = 0.7, and 
H = 70 km s _1 Mpc -1 . 



2 THE COSMIC HORSESHOE 

The Cosmic Horseshoe, or Jl 148+1930, is a nearly complete 
blue Einstein ring of 10 arcsec diameter encircling a red galaxy 
(FigureQ] see also http://www.ast.cam.ac.uk/research/cassowary/). 
Low resolution spectroscopy by Belokurov et al. (2007) revealed 
the lens to be a massive luminous red galaxy (M ~ 6 x 1O 12 M0) 
at z — 0.444 and the ring to be the gravitationally lensed image of 
a star-forming galaxy at z = 2.379. From the extent of the ring, 
Belokurov et al. estimated an approximate magnification factor of 
~ 35, while the detailed lensing models constructed by Dye et al. 
(2008) give overall magnification factors of 24 ± 2. Adopting the 
latter value, and interpolating between the total integrated magni- 
tudes of the ring g = 20.1, i = 19.7 reported by Belokurov et al. 
(2007), we deduce a luminosity L ~ 2.4L* for the source, rela- 
tive to MX B (1700 A) = -21.0 from the luminosity function of 
UV-selected z ~ 2 galaxies by Reddy et al. (2008). The image re- 
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construction in the source plane with the best fitting lensing model 
by Dye et al. (2008) shows that the major contributor to this lumi- 
nosity is a compact source with half-light radius rj/2 ~ 0.3 arcsec, 
or ~ 2.5 kpc at z = 2.379; there is also evidence for a second, 
fainter, source ~ 0.7 arcsec to the North. 

We briefly review the information currently available on the 
lensed galaxy. Its rest-frame UV spectrum is essentially flat in F v , 
reflecting a recent star-formation episode and modest reddening 
[E(B - V) ~ 0.1-0.2; Belokurov et al. 2007], but a more ex- 
tensive characterization of its stellar population(s) awaits longer 
wavelength photometry to build-up the spectral energy distribution 
from the UV to the near-IR (e.g. Shapley et al. 2005; Erb et al. 
2006b). Recently, Hainline et al. (2009) reported observations of 
several rest-frame optical emission lines formed in the H II regions 
of the Cosmic Horseshoe, from which they were able to measure 
a number of parameters of interest. The Ha luminosity implies a 
star formation rate of 113 M@ yr _1 , after correcting for reddening 
and a lensing magnification factor of 24. The emission line widths 
indicate a velocity dispersion a ~ 65 km s" 1 and a correspond- 
ing virial mass M v j r ~ 1.0 x 1O 1O M0. The oxygen abundance 
deduced from the ratios of the strongest emission lines is in the 
range (0/H)h ii — 0.5-1.5 (O/H)©. These parameters are typical of 
UV-bright galaxies at z ~ 2 (Erb et al. 2006a,b,c), although the 
star formation rate is at the upper end of the distribution of values 
found by Erb et al. (2006c). It is of interest now to compare such 
estimates with those provided by the rest-frame UV spectrum. 



3 OBSERVATIONS AND DATA REDUCTION 

We observed the Cosmic Horseshoe with the Echellette Spectro- 
graph and Imager (ESI; Sheinis et al. 2000) on the Keck II tele- 
scope. The high efficiency, moderately high resolution (R — 4000 
with a 1 arcsec slit, sampled with 11.4 km s" 1 pixels), and wide 
spectral coverage of ESI— from 4000 to 10000 A, or 1184-2959 A 
in the rest-frame of the Cosmic Horseshoe — make it highly suit- 
able for an in-depth study of the UV spectrum of this galaxy. The 
data were collected over two nights (2008 March 6 and 7 UT) in 
mostly sub-arcsecond seeing; the total exposure time was 36100 s, 
made up of a number of 1800 s or 2000s-long individual integra- 
tions. The ESI slit was oriented at sky RA. = 10.8 degrees East of 
North, and positioned so as to encompass the two brightest knots in 
the Einstein ring, as shown in Figure[T] The northern knot, labelled 
'Aperture 1 ' in the Figure, is knot 'A' in the labelling of the features 
by Belokurov et al. (2007), while Aperture 2' corresponds to their 
knot 'D'. Both knots are due to the same object in the source plane, 
although Aperture 1' may include a small contribution from the 
fainter component ~ 0.7 arcsec to the North of the main source of 
the ring (Dye et al. 2008). The location of the ESI slit is the same 
as that of the NIRSPEC observations by Hainline et al. (2009). 

The data were processed with standard IRAF tasks. Each in- 
dividual two-dimensional (2-D) ESI spectrum was geometrically 
corrected (to correct for the spectral curvature of each echellette 
order), bias subtracted, flat-fielded, and corrected for cosmetic de- 
fects and cosmic-ray affected pixels. In the next step, all of the 
images thus processed were co-added into a total 2-D frame which 
was then background subtracted. The background subtraction pro- 
cess included careful modelling of the distribution along the slit of 
light from the lens, which is clearly visible above the sky at the 
longer wavelengths of the ESI range. Finally, two 1-D spectra were 
extracted (using optimal extraction algorithms), one for each 'Aper- 
ture'. The extraction process included merging the 10 echellette or- 



ders of ESI into one continuous spectrum, with weighted addition 
of portions of adjacent orders which overlap in wavelength. 

We compared carefully the 1-D spectra of Aperture 1 and 
Aperture 2 to look for differences which one may expect to be re- 
vealed by the gravitational magnification of any spatial structure 
in the source, but found none. Neither dividing one spectrum by 
the other nor subtracting one from the other revealed features more 
significant than the noise. As discussed in Section[6] this is also the 
case for the Lya emission line which (a) has the highest signal-to- 
noise ratio (S/N), and (b) is most sensitive to geometrical effects, 
being resonantly scattered. The only difference between the two 
apertures is in the total flux: emission line and continuum are uni- 
formly higher by a factor of 1.19 in Aperture 2 than in Aperture 1. 
Thus, in order to improve the S/N, the two spectra were averaged 
and the resulting spectrum mapped onto 0.5 A bins. 

This final spectrum has a resolution (full width at half maxi- 
mum) FWHM = 75 km s ~~ 1 , sampled with three wavelength bins at 
6000 A, as determined from the widths of narrow emission lines 
from the Cu-Ar and Hg-Ne-Xe hollow-cathode lamps used for 
wavelength calibration, whose spectrum was processed in the same 
way as that of the Horseshoe. From the rms deviation of the data 
from the mean, we measured an average S/N ~ 8 per 0.5 A wave- 
length bin between ~ 4000 and ~ 7500 A (~ 1200 and ~ 2200 A 
in the rest frame), with a factor of ~ 50% variation in S/N across 
each echellette order. At the wavelengths of the Lya emission line 
the S/N rises to a maximum S/N = 44 per 0.5 A bin. 

We attempted to put our spectrum on an absolute flux scale 
by reference to those of three flux-standard stars recorded dur- 
ing the two nights of observation, but found a scatter of ±20% 
in the flux calibration among the three stars. Adopting the mean 
value of the three flux scales, we deduced a mean (/„) = 
8.3 x 10~ 29 erg s _1 cm" 2 Hz" 1 between 4200-5485 A. This 
wavelength interval corresponds to the FWHM of the transmis- 
sion curve of the Gunn g filter through which Belokurov et al. 
(2007) measured g = 20.1 (corresponding to {/„) = 3.3 x 
10 -28 erg s _1 cm -2 Hz -1 ) for the entire Einstein ring. Thus, 
if the flux calibration is correct, the ESI slit captured a fraction 
0.25 ± 0.05 of the light of the Cosmic Horseshoe, corresponding to 
an effective magnitude within the two Apertures of g = 21. 6± 0.2. 



4 THE STELLAR SPECTRUM 

The rest-frame UV spectra of star-forming galaxies are complex 
blends of interstellar absorption lines, nebular emission lines, and 
absorption/emission lines formed in the atmospheres of OB stars. 
This wealth of spectral features, properly interpreted, is a rich 
source of information on the physical properties of the gas and stars 
in these galaxies. In this section, we consider photospheric absorp- 
tion lines which have been shown to be useful abundance diagnos- 
tics, and P-Cygni profiles formed in the expanding atmospheres of 
the most luminous early-type stars whose strength is also sensitive 
to the upper end of the initial mass function. 

4.1 The systemic redshift of the Cosmic Horseshoe 

The first step in the analysis is the determination of the systemic 
redshift of the galaxy. While most photospheric features are un- 
resolved blends of multiple lines, a few exceptions allow the red- 
shift of the hot stars to be determined. Within the wavelength range 
covered by our data, we identified the following unblended pho- 
tospheric absorption lines: SVA1501.763, N IV A1718.551, and 
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Table 1. Systemic Redshift 



Ion 


A fab <A) 


-^•abs 


Origin 


S V 


1501.763 


2.38107 


Stars 


NIV 


1718.551 


2.38119 


Stars 


Cm] 


1906.683 


2.38103 


H II Regions 


Cm] 


1908.734 


2.38127 


H II Regions 


C III 


2297.579 


2.38118 


Stars 


Ha b 


6564.614 


2.38123 


H II Regions 



a Vacuum wavelengths 

b As reported by Hainline et al. (2009) 



CIIIA2297.579 [A air = 2296.871— all three vacuum wavelengths 
from the National Institute of Standards and Technology (NIST) 
database available at http://physics.nist.gOv/PhysRefData/ASD/l. 

The redshifts of these three lines are in good mutual agreement 
(see Table[T), and define an average z s tars = 2.38115 ± 0.00006. 
The la error corresponds to an uncertainty of ±5km s _1 , or 1/5 
of a wavelength bin in our final spectrum. Also included in TablefT] 
are values of z Hu from the C III] AA1906.683, 1908.734 doublet 
lines which are resolved in our ESI spectrum (see Section [7J, and 
from Ha measured from the NIRSPEC spectrum of Hainline et 
al. (2009). The mean z Cm] = 2.38115 ± 0.00012 is in excellent 
agreement with z B t alB ', zu a differs by only 7km s _1 which may 
include a small systematic offset between the wavelength calibra- 
tions of the ESI and NIRSPEC spectra. Thus, we conclude that the 
systemic redshift of the Cosmic Horseshoe is jz sys = 2.38115 and 
we adopt this value throughout the paper. 

4.2 Photospheric lines 

In this and the following subsection, we analyse photospheric and 
wind lines in the spectrum of the Cosmic Horseshoe by compar- 
ing our data with model spectra computed with the population syn- 
thesis code Starburst99 which couples libraries of either empirical 
(Leitherer et al. 1999, 2001) or theoretical (Rix et al. 2004) ultravio- 
let OB stellar spectra with stellar evolutionary tracks. The input set 
of parameters to Starburst99 can be adjusted to simulate a variety 
of star formation histories, ages, IMF parameters, and metallicities 
so as to determine the combination that best fits the observations 
under scrutiny. 

In the absence of more detailed information on the previous 
history of star formation in the Cosmic Horseshoe, we shall adopt 
the simplest case of continuous star formation — with a Salpeter 
slope of the upper end of the IMF — which has been proceeding 
at a steady rate (as opposed to a declining rate, or a series of bursts) 
for 100 Myr (this being a safe lower limit beyond which the spec- 
trum no longer changes with time because the stars that contribute 
to the UV light are born and die at the same rate). While this sce- 
nario is undoubtedly an over-simplification, it may not be too differ- 
ent from reality when considering the spectrum of a whole galaxy, 
rather than individual regions of star formation. An age of more 
than 100 Myr is typical of most z — 2 UV-bright galaxies in the 
'BX' sample of Erb et al. (2006b) and seems plausible given the 
near-solar metallicity of the H II regions of the Horseshoe (Hain- 
line et al. 2009). 

Leitherer et al. (2001) and Rix et al. (2004) pointed out the ex- 
istence of some blends of UV photospheric lines whose strengths, 
under the above assumptions, depend on metallicity (see also Hal- 
liday et al. 2008). The "1425" index of Rix et al. (2004) mea- 



sures the equivalent width of a blend of Si III A1417, Cm A1427, 
and Fe V A 1430 absorption lines spanning the wavelength interval 
1415-1435 A; at longer wavelengths, the "1978" index measures 
the strengths of several Fe III absorption lines from B stars between 
1935 and 2020 A. Both indices stabilise after ~ 50 Myr from the 
onset of star formation, and increase monotonically with metallic- 
ity. In the left-hand portion of Figure [2] we compare the ESI spec- 
trum of the Cosmic Horseshoe with the model spectra of Rix et 
al. (2004) for five values of metallicity, from 1/20 of solar to twice 
solar, in the regions of the 1425 and 1978 indices. For this compari- 
son, we smoothed our spectrum to the 2.5 A resolution of the Rix et 
al. models, and normalized it using the psuedo-continuum windows 
suggested by those authors. Thus, the observed and model spectra 
should be directly comparable. 

Figure|2]shows a generally good agreement in the detailed fea- 
tures of the synthetic and real spectra in the two wavelength regions 
considered here, attesting to the sophistication reached by current 
hot-star atmosphere and stellar population synthesis models. Re- 
ferring to the 1425 A region (left-most panel in Figure|2]l, it can be 
seen that the strength of the photospheric lines in this region of the 
Cosmic Horseshoe spectrum is intermediate between those com- 
puted with metallicities Z = QAZq and Z — I.OZq respectively. 
Using the relationship between metallicity and the 1425 index pro- 
posed by Rix et al. (2004): 

log(Z/Z Q ) —Ax EW(1425) + B (1) 

with A = 1.14 and B — —1.75 appropriate to the value 
EW(1425) = 0.95 A we measure here, we derive ZoBstars = 
0.5Z Q . 

Qualitatively, the observed spectrum seems to be intermedi- 
ate between the 0.4 solar and solar metallicity cases in the 1978 
region too (see middle panel of Figure[2]l. However, there appears 
to be excess absorption between 1960 and 1980 A over that ex- 
pected from the models for any metallicity; the presence of these 
additional features does not allow us to deduce a value of Zestars 
from EW(1978) using the Rix et al. (2004) calibration of this index. 
We considered the possibilities that the excess absorption may be 
noise, or an intervening absorption system unrelated to the Cosmic 
Horseshoe, but could find no evidence to support either interpre- 
tation. It also seems unlikely that the excess is due to a different 
mix of stellar spectral types or chemical elements, since all of the 
photospheric absorption in this region is thought to be due to Fe III 
lines, and the rest of the spectral features between 1920 and 2025 A 
match the data well. No such discrepancies were found by Rix et 
al. (2004) in the two galaxies they considered (MS 1512-cB58 and 
Q1307-BM1163, the latter at z = 1.411), nor by Halliday et al. 
(2008) in their composite spectrum of z ~ 2 galaxies. However, 
the cases where these photospheric indices have been measured are 
still very few, and it is important to continue to look critically at 
the match between model and real spectra as more data of suitable 
quality become available. 



4.3 Wind lines 

The model spectra by Rix et al. (2004) show that the spec- 
tral features most sensitive to metallicity are the P-Cygni lines 
formed in the expanding winds of the most luminous OB stars 
whose mass-loss rates are thought to be lower at lower metal- 
licities (e.g. Kudritzki & Puis 2000). Among the P-Cygni lines 
covered by our spectrum, the C IV AA1548.204, 1550.781 (un- 
resolved) doublet is the strongest. Referring to the rightmost 
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Figure 2. Left-hand pair of panels: Portions of the UV spectrum of the Cosmic Horseshoe (black) in the regions encompassing the 1425 and 1978 indices 
compared with model spectra (red) from the work by Rix et al. (2004) for five different metallicities, as indicated. The ESI spectrum has been smoothed to the 
2.5 A resolution of the model spectra. See text for further details of the models. Right-hand panel: Comparison of the Rix et al. models with the smoothed and 
normalized spectrum of the Cosmic Horseshoe in the region of the C IV AA1548.204, 1550.781 doublet. This region also includes C IV and Si II Af 526.7070 
interstellar absorption which is unaccounted for in the Rix et al. models. Note the difference in the j/-axis scales between the C IV panels and those for the 
weaker photospheric 1425 and 1978 features to the left. 



panel in Figured we find that its strength in the Cosmic Horse- 
shoe is intermediate between those of the models with metal- 
licities Z = 0.2Zq and Z = OAZq. However, as empha- 
sised by Crowther et al. (2006), the interpretation of this feature 
is complicated by its blending with interstellar C IV (and to a 
lesser extent Si II A1526.7070) absorption, at least at the coarse 
resolution of most spectra of distant galaxies. Similar problems 
affect the generally weaker Si IV AA1393.7602, 1402.7729 and 
N V A1238.821, 1242.804 wind lines. 

At the high resolution of our ESI spectrum, interstellar and 
stellar components are easily separated, as can be appreciated from 
Figure [3] The top panel compares the observed spectrum of the 
Cosmic Horseshoe in the C IV region with those generated by Star- 
burst99 now using libraries of empirical stellar spectra, assem- 
bled from UV observations of stars in either the Galaxy or the 
Magellanic Clouds. Both models are for 100 Myr old continuous 
star formation with Salpeter IMF. Apart from the narrow interstel- 
lar SillA1526.7070 and C IV AA1548.204, 1550.781 absorption 
lines, which are stronger in the distant galaxy than in the spec- 
tra of the stars that make up the Starburst99 libraries, the agree- 



ment between models and observations is clearly very good. The 
broad blue absorption wing of the P-Cygni profile reaches simi- 
lar terminal velocities in the models and the data, and its depth is 
roughly intermediate between those of the Galactic and Magellanic 
Clouds models. The latter library of stellar spectra was built-up by 
Leitherer et al. (2001) with Hubble Space Telescope (HST) obser- 
vations of stars in both the Large and Small Magellanic Clouds 
(LMC/SMC) and may thus correspond to an approximate metal- 
licity Zmc ~ 0.4Z Q , given that in the LMC (0/H)~ O.5(O/H) 
and in the SMC (O/H) ~ 0.25 (O/H) (e.g Pagel 2003). By lin- 
early interpolating between the Galactic and LMC models and 
minimising the difference from the observed spectrum, we find 
Zostars « 0.6Z Q , in good agreement with ZoBstars = 0-5Z e 
deduced above from consideration of the 1425 photospheric index. 

It is interesting that both Galactic and Magellanic Clouds 
models seem to underpredict the emission component of the P- 
Cygni profile. However, the narrow emission feature centred near 

155 1.6 A may be partly of nebular origin (Leitherer, Calzetti, 
& Martins 2002). Similarly, the narrow emission centred near 

1533.7 A is probably nebular Sill* A1533.4312; both features are 
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Figure 3. Upper panel: Comparison between the ESI spectrum of the Cosmic Horseshoe in the region encompassing the C IV line and model spectra computed 
with Starburst99 and empirical libraries of Galactic and Magellanic Clouds stars, as indicated. The model spectra were generated assuming a 100 Myr old 
continuous star formation episode with a Salpeter IMF. Lower panel: The stellar spectra of the Cosmic Horseshoe and of MS 1512-cB58 are remarkably 
similar in the wavelength region shown. The ESI spectrum of cB58 is reproduced from Pettini et al. (2002), and has been reduced to its rest wavelengths at 
Zstars = 2.7276. 



missing from the synthetic spectra which were designed to repro- 
duce only the stellar component of the galaxy spectrum. The super- 
position of P-Cygni broad emission/absorption, photospheric broad 
absorption, narrow interstellar absorption and narrow nebular emis- 
sion attests to the complex nature of this portion of the UV spec- 
trum of star-forming galaxies. 

In the lower panel of Figure [3] we compare the ESI spectra of 
MS 1512-cB58 (from Pettini et al. 2002) and the Cosmic Horseshoe 
in the wavelength region 1500-1570 A. There is a startling similar- 
ity in this spectral region between these two galaxies which are at 
different redshifts (cB58 is at z s tars = 2.7276) and were selected 
randomly, only by virtue of the fact that they are highly lensed as 
seen from Earth (although they are both more luminous than the 
'average' galaxy at these redshifts, even after correcting for the 
lensing magnifications). Presumably, not only the chemical abun- 
dances, but also the young stellar populations of these two galaxies 
are remarkably similar. Only the interstellar absorption lines are 
different, with those in cB58 apparently stronger than in the Cos- 
mic Horseshoe (but see the discussion in Section[5]l- 

The P-Cygni component of the C IV complex is due to the 
most massive (and luminous) O stars, which drive the strongest 
winds. On the other hand, a wider range of stellar spectral types, in- 
cluding all stars with masses greater than about 5A-/0 , are thought 
to make up the integrated continuum light near 1550 A (e.g. Rix et 



al. 2004). Thus, the contrast of the P-Cygni features relative to the 
continuum is sensitive not only to metallicity, but also to the slope 
and upper end cut-off of the IMF, as well as to more subtle effects, 
such as differential dust extinction among stars at the upper end of 
the IMF (the most massive and short-lived stars may in principle 
be more reddened than longer lived ones — Leitherer et al. 2002), 
and the relative proportions of OB stars in star clusters and in the 
field (Chandar et al. 2005). Pettini et al. (2000) showed how even 
relatively minor changes to the IMF slope (or alternatively the up- 
per mass cut-off) result in noticeable alterations to the integrated 
spectrum of a star-forming galaxy in the C IV region. On that ba- 
sis, those authors concluded that a standard Salpeter IMF extend- 
ing to A/ U p > 50Mq provides the best match to the spectrum of 
cB58, without the need to invoke a different IMF at high redshift, 
as claimed by some. Clearly, the same conclusion applies here to 
the massive stellar population in the Cosmic Horseshoe. More gen- 
erally, the similarity evident in the lower panel of Figure[3]suggests 
that the minor effects discussed above (differential dust extinction 
and the balance between cluster and field stars) presumably aver- 
age out in at least a subset of star-forming galaxies at z = 2 — 3, 
when their integrated spectra are considered. 

Summarising the conclusions of this section, we find that 
the metallicity of the early-type stars in the Cosmic Horseshoe is 
approximately half-solar, from consideration of photospheric and 
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Figure 4. Normalized profiles of selected interstellar absorption lines. Left panel: Transitions of ions which are the dominant species of the corresponding 
elements in H I regions. Right panel: Higher ionization stages. In both panels the vertical red dash lines indicate the velocity range over which the values of 
equivalent width listed in Table|2]were generally measured. 



wind absorption lines. A continuous mode of star formation with a 
Salpeter slope for stars more massive than 5 Mq gives a good rep- 
resentation of the UV spectrum. The 1978 index introduced by Rix 
et al. (2004) appears to be contaminated by additional features of 
uncertain origin in the spectrum of the Horseshoe. The C IV region 
is a complex blend of stellar, interstellar and nebular features, but 
the stellar component in the Cosmic Horseshoe is surprisingly sim- 
ilar to that in cB58, the only other high redshift star-forming galaxy 
studied in comparable detail so far. This may of course be a coinci- 
dence, but it may also point to the fact that many of the parameters 
which can affect the appearance of this spectral region average out 
to a large extent in the integrated UV spectrum of a whole galaxy. 



5 THE INTERSTELLAR SPECTRUM 
5.1 Kinematics of the absorbing gas 

In the wavelength interval from 1200 to 2600 A, we identified 21 
interstellar absorption lines from six chemical elements in a variety 
of ionization stages, from O I to C IV. Relevant measurements are 
collected in Table[2] while Figure|4]shows selected transitions cho- 
sen to illustrate the range of elements and ion stages covered. The 
interstellar lines in the Cosmic Horseshoe are broad, with absorp- 
tion extending smoothly from —800 to +250 km s _1 with respect 
to z sya = Zstars = 2.38115. Within the limits of the noise, we find 
no convincing differences between the profiles of different ions, al- 



though the maximum velocities measured depend on the strength 
of the transition, as is the case in local starbursts (e.g. Grimes et al. 
2009). The values of z a bs listed in column (4) of Table[2]are those 
derived from the centroids of the lines (that is, the mean wavelength 
of the line weighted by the absorption in each wavelength bin); their 
mean value is (z) = 2.3788 ± 0.0009 which corresponds to a ve- 
locity offset Av = —208 km s -1 with respect to the stars. We also 
measured the redshift at the peak optical depth in the most clearly 
defined line profiles and found a mean zism = 2.3795 ± 0.0002 
(Av = -146 kms" 1 ). 

This net blueshift of the interstellar lines is a common feature 
of star-forming galaxies at low (e.g. Heckman et al. 2000; Martin 
2005) as well as high (e.g. Pettini et al. 2001; Shapley et al. 2003; 
Steidel et al. 2009; Vanzella et al. 2009) redshifts. It is generally 
interpreted as empirical evidence for the existence of large-scale 
outflows of the interstellar medium driven by the kinetic energy 
deposited by supernovae and the winds of massive stars. The values 
of Av we deduce here, irrespective of whether one considers the 
mean redshift of the line centroids or of the gas with the highest 
optical depth, are typical of z — 3 galaxies for which Shapley et al. 
(2003) derive a mean (Aw) = -150 ± 60km s~\ Similarly, the 
widths of the strongest lines (FWHM ~ 400-600 km s _1 ) are 
within the range (FWHM) = 575 ± 150 km s _1 measured by 
Shapley et al. (2003) from their composite spectrum of 81 1 LBGs. 

Although with the higher resolution of our data (compared to 
that normally employed to record the spectra of unlensed z — 2 — 3 
galaxies) the interstellar lines are fully resolved, we are still unable 
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Table 2. Interstellar Absorption Lines 



Ion 


Ai ab a (A) 


/ a 




W h (A) 


<x b (A) 


Comments 


Cn 


1334.5323 


0.1278 


2.3786 


1.62 


0.16 


Blended with C II* A1335.6627 


Civ 


1548.204 


0.1899 


2.3783 


2.36 


0.08 


Blended with stellar C IV A1549.1; Wo measured from 




1550.781 


0.09475 


2.3798 


0.85 


0.04 


Blended with stellar C IV A1549.1; Wo measured from 


Ol 


1302.1685 


0.04887 


2.3791 


0.90 


0.09 


Wo measured from —800 to +100 km s -1 


Alii 


1670.7886 


1.74 


2.3779 


1.48 


0.07 




Aim 


1854.7184 


0.559 


2.3792 


0.93 


0.10 






1862.7910 


0.278 


2.3782 


0.61 


0.13 




Sin 


1260.4221 


1.18 


2.3778 


1.71 


0.18 


Blended with S II A1259.519 




1304.3702 


0.0863 


2.3792 


0.71 


0.06 


Wo measured from —400 to +25 km s — 1 




1526.7070 


0.133 


2.3788 


1.81 


0.06 






1808.0129 


0.00208 


2.3778 


0.40 


0.08 




Si IV 


1393.7602 


0.513 


2.3795 


1.36 


0.09 






1402.7729 


0.254 


2.3791 


1.13 


0.08 




Fen 


1608.4511 


0.0577 


2.3774 


0.51 


0.12 






2344.2139 


0.114 


2.3807 


0.40 


0.23 






2374.4612 


0.0313 


2.3795 


1.47 


0.23 






2382.7652 


0.320 


2.3767 


1.25 


0.18 


Wo measured from —700 to +250 km s _1 




2586.6500 


0.0691 


2.3791 


1.85 


0.13 






2600.1729 


0.239 


2.3793 


2.25 


0.28 




Mil 


1317.217 


0.0571 


2.3786 


0.12 


0.05 


Wo measured from —350 to —80km s — 1 




1741.5531 


0.0427 


2.3799 


0.14 


0.05 


Wo measured from —180 to —20km s _1 



-800 to +100kms _1 
-250 to +60kms" 1 



a Vacuum wavelength and /-values are from Morton (2003) with updates by Jenkins & Tripp (2006). 

b Rest-frame equivalent width and associated error measured over the velocity range —800 to +250 km s _1 with respect to 2 s tars, unless otherwise 
noted. 



to recognize clear clues in the line profiles as to the location of the 
absorbing gas. The profiles are smooth and relatively featureless, 
and we see no trends between, for example, velocity and degree of 
ionization, which could be used to construct a kinematical model. 
While the blueshifted absorption is presumably due to outflowing 
material, the location and nature of the gas moving at positive ve- 
locities relative to the stars remain unexplained. It is interesting to 
note that the absorption lines in the ESI spectrum of MS 1512-cB58 
analysed by Pettini et al. (2002) also extend over the same velocity 
range, from ~ —800 to ~ +250 km s _1 relative to the systemic 
redshift of the stars, although in cB58 the gas with the highest opti- 
cal depth is moving at a higher velocity Av — —255 km s . The 
two galaxies have similar star formation rates, SFR ~ 50 Mq yr~ 1 , 
as estimated from their UV continua (see Section [o3t , Overall, the 
profiles of the interstellar absorption lines from the first ions appear 
smoother in the Horseshoe than in cB58, but the lower S/N of the 
present data makes it more difficult to recognize distinct velocity 
components than is the case in cB58. 



5.2 Partial coverage and column densities 

One noticeable difference between the Cosmic Horseshoe and 
cB58 is in the optical depth of the interstellar absorption lines (see 
Figure(3](- Whereas in cB58 the strongest interstellar lines have sat- 
urated cores of zero residual intensity (at the spectral resolution of 
ESI — see Figures 1 and 2 of Pettini et al. 2002), the same tran- 
sitions never seem to reach below an optical depth r ~ 1 in the 
Cosmic Horseshoe, as can be appreciated from inspection of Fig- 
ure [4] 

One possibility is that the column densities of even the most 
abundant ions are genuinely lower in the Horseshoe compared to 
cB58. An alternative explanation is that the interstellar gas does not 
completely cover the early-type stars producing the UV continuum 



against which the absorption is seen. The latter interpretation is not 
unlikely in general — given the composite nature of the UV spectra 
of star-forming galaxies which are the superposition of hundreds 
of thousands of individual stellar spectra — and is suggested by the 
similarity in the residual intensity in the cores of all the strongest 
absorption lines in Figure|4] 

We can assess quantitatively the partial coverage hypothesis 
using the apparent optical depth method of Savage & Sembach 
(1991). If the absorption lines are resolved, as seems to be the case 
here, the column density of an ion in each velocity bin, N a (v) [in 
units cm -2 (km s -1 ) -1 ], can be deduced directly from the optical 
depth in each velocity bin, 

T«,(«)=-In[J b.(t»)//o(«)] (2) 

from the relation 

N a (v) = 3.768 xlO 14 -^^ (3) 

J A 

where / measures the strength of the transition at wavelength A 
(A), and / b s and Io denote the relative intensities in the line and 
in the continuum respectively. If multiple absorption lines arising 
from the same ground state of an ion but with different values of the 
product /A are available, partial coverage would manifest itself as a 
mismatch between the values of column density deduced from each 
transition, the lines with smaller values of /A giving systematically 
higher values of N a (v). 

Among the ions covered by the ESI spectrum of the Horse- 
shoe, Si II is the one best suited to this analysis, with four transi- 
tions spanning a range of ~ 400 in /A (see Table |2) and falling at 
wavelengths where the S/N of the data is highest. As can be seen 
from Figure [5] the run of N a (v) vs. v for the four Sill lines does 
indeed show the systematic differences expected if the absorbing 
gas does not completely cover the background stellar continuum. 
In the line core, the apparent optical depth of the weakest absorp- 
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Figure 5. When analysed with the Apparent Optical Depth method, the four 
Si II transitions recorded in our ESI spectrum of the Cosmic Horseshoe yield 
values of column density which differ systematically, in the sense expected 
if the interstellar gas does not completely cover the stellar UV continuum 
(see Section l531 for further details). 



tion line, Si II A1808.0129, implies values of column density two 
orders of magnitude higher than those deduced from the strongest 
transition, Si II A1260.4221. The profiles of Si II A1304.3702 and 
A1526.7070, with very similar values of /A intermediate between 
those of the other two lines, indicate that partial coverage extends 
well beyond the line cores and probably applies to the full velocity 
range spanned by the absorption. 

We also cover several Fell transitions (see Table[2j, but their 
/A values range over a factor of only ~ 10 and most of them 
are redshifted to far-red wavelengths, where the S/N of our data 
is lower. Nevertheless, the profiles of the Fell lines are consis- 
tent with the conclusions drawn from Figure [5] in showing a sys- 
tematic trend of increasing N a (v) with decreasing /A. Turning to 
the high ionization lines in the right-hand panel of Figure [4] we 
found that the weaker member of the Si IV doublet, A1402.7729, 
also gives systematically higher values of N a (v) than A1393.7602, 
even though the /A values of the two lines differ by only a factor 
of 2, while the evidence is less clear-cut in the case of the Al III 
doublet which is recorded at lower S/N. The C IV lines are difficult 
to interpret in this context because they are blended with each other 
and with the stellar P-Cygni profile, as discussed above (Section 

Given our lack of knowledge of the relative configuration of 
interstellar gas and stars, it is quite possible, or even likely (e.g. 
Martin & Bouche 2009), that the covering factor varies with veloc- 
ity and degree of ionization of the gas. However, the S/N of our data 
is insufficient to establish the magnitude of such variations. Instead, 
we adopted the simplest assumption that "one size fits all" and cor- 
rected the absorption line profiles by re-adjusting the zero level by 
the same fractional amount of the continuum at all wavelengths. In- 
spection of the spectrum shows that even the strongest absorption 
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Figure 6. Black histograms: Normalized profiles of the four Si II absorp- 
tion lines after subtracting 40% of the continuum level, under the assump- 
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and velocity dispersion parameter b implied by the weakest transition, 
A1808.0129, which is unsaturated. The absorption line shown by the short- 
dash histogram in the top right-hand panel is O I A1302. 1685. 



lines, which are normally saturated in the spectra of Galactic stars 
and of nearby starbursts, all seem to have a minimum residual in- 
tensity Ix/Io ~ 0.4 (see FigureHJ. With the simplistic assumption 
that the interstellar gas only covers 60% of the stellar light, and 
that the remaining 40% of the UV continuum sees no absorption 
in our direction, the adjusted profiles of the four Si II lines become 
internally consistent, as shown in Figure[6] 

Specifically, after subtracting 40% of the stellar continuum 
and renormalizing the line profiles, we deduced the column den- 
sity of Si II, log iV(Si Il)/cm~ 2 = 16.00, by integrating eq. l|3j over 
the velocity range spanned by the weakest Sill line, A1808.0129, 
which is unsaturated even after the readjustment of the zero level. 
We then used the absorption line fitting software XVOIGT (Mar & 
Bailey 1995) to fit the profile of Si II A1808.0129 with this column 
density and deduced the value of the velocity dispersion parameter, 
b = 105 km s _1 , which best reproduces the width of this line. In 
the next step, we used this pair of values of iV(Si II) and b to gen- 
erate XVOIGT profiles of the other three Si II absorption lines. As 
can be seen fom Figure [6] these computed profiles provide satis- 
factory fits to the main absorption component of the three stronger 
Si II transitions, given the noise in the data. (We have not attempted 
to fit the blue wing of the three lines because this absorption is too 
weak to be detected in Si II A1808.0129, but its contribution to the 
value of A f TOT(Si II) is only a small fraction of the total, unless the 
covering factor of the high velocity gas is much smaller than 60%). 

Table [3] lists the values of ion column density deduced by 
integrating eq. over the velocities spanned by those absorp- 
tion lines that remain unsaturated after the 40% adjustment of the 
zero level. Errors in the column density, 8N, were determined by 
adding/sub trading the la error spectrum to/from the line profiles 
and recalculating the column density; in cases where the addition 
of the noise spectrum results in the line becoming saturated no cor- 
responding error is given in Table [3] Again, we see from the Table 
that there is good internal agreement between the values of N im- 
plied by different transitions of the same ion, after correcting for 
the 40% unabsorbed stellar continuum. 
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Table 3. Ion Column Densities 



Ion 


Alab (A) 


/ 


log Af/cm 2 


log 5/V/cm 2 


Aim 


1854.7184 


0.559 


14.22 


-0.06, ... 




1862.7910 


0.278 


14.26 


-0.07, +0.14 


Sin 


1808.0129 


0.00208 


16.00 


-0.06, +0.07 


Si IV 


1402.7729 


0.254 


14.92 


-0.06, ... 


Fell 


1608.4511 


0.0577 


14.94 


-0.10, +0.11 




2586.6500 


0.0691 


15.12 


-0.06, ... 


Nill 


1317.217 


0.0571 


14.22 


-0.24, +0.20 




1741.5531 


0.0427 


14.40 


-0.21, +0.19 



The interstellar column densities of the first ions in the 
Cosmic Horseshoe are surprisingly similar to those measured 
in cB58 by Pettini et al. (2002). For the transitions in com- 
mon between the two studies, we have: log TV/cm -2 (Horse- 
shoe/cB58) = 16.00/15.99 for Si II A1808.0129; 14.94/15.17 
for FenA1608.4511 ; 14.22/14.32 for Ni II A1317.217; and 
14.40/ < 14.45 for Ni II A1741.5531. Although the metallicities 
reached by the two galaxies are similar, as discussed above (Sec- 
tion[4]l, there was no reason to expect that the column densities of 
interstellar gas entrained by the outflows — nor the ionization of this 
gas — should be so similar. Furthermore, even the ratios of elements 
that are depleted onto dust by different degrees, or are synthesised 
by stars of different masses, are similar between these two galax- 
ies. Thus, the iron-peak and easily depleted elements Fe and Ni are 
underabundant (in their first ions) compared with the undepleted 
alpha-capture element Si by about one order of magnitude, as is the 
case in cB58 (see Figure 8 of Pettini et al. 2002). All of this may 
just be a coincidence, of course, or may point to a more fundamen- 
tal uniformity in the conditions of the gas outflowing from high- 
redshift star-forming galaxies, or at least a subset of such galaxies. 

Given the similarity in the column densities of the first ions 
and in the metallicity of the two galaxies, it follows that we also 
expect similar values of the column density of hydrogen 7V(H). In- 
deed, if the majority of the first ions arises from Hi regions, we 
can deduce values of iV(H I) from the column densities in Table|3] 
and the assumption that Zism = ZoBstars = 0.5Zq. This leads to 
log 7Y(H l)/cm -2 = 20.79 usina Si II as the reference ion, indeed 
very close to log 7V(H l)/cm~ 2 = 20.85 deduced by Pettini et al. 
(2002) from the damped profile of the Lya line in cB58. As we shall 
see in a moment, however, there are striking differences between 
the profiles of the Lya line observed in the two galaxies. Whereas 
in the case of cB58 we see a clear damped absorption profile with 
a weak superposed emission component, Lya is in emission in the 
Cosmic Horseshoe. 



6 THE LYMAN ALPHA LINE 
6.1 Lyman a line morphology 

The top panel of Figure [7] shows the spectral region encompassing 
the Lya line in the Cosmic Horseshoe, plotted in arbitrary flux- 
density units vs. velocity relative to the systemic redshift of the 
galaxy. We interpret the complex morphology as consisting of two 
main emission components: a narrow, symmetric component which 
is barely resolved (FWHM 50 km s _1 after deconvolution with 
the ESI spectral resolution of FWHM = 75 km s ) centred at v — 
+ 115 km s _1 , and a broader, asymmetric component which peaks 
at v = +275 km s" 1 but extends to v ~ +700 km s _1 . The first, 



narrow component has a peak flux ~ 1.5 times higher than that 
of the second. There may be a third, weak component at negative 
velocities, centred near v ~ —80 km s _1 , but its identification is 
uncertain, given the blending with absorption from the Lya forest. 
The total flux in the line, without correcting for gravitational lens- 
ing and reddening (these corrections will be applied later — see Sec- 
tionEOt is F(Lya) = (4.5+0.2) x 10" 16 erg s" 1 cm" 2 where the 
error reflects both the random noise in the data and the uncertainty 
in the stellar continuum (measured redward of Lya), combined in 
quadrature. In our assumed cosmology, the measured flux corre- 
sponds to aLya luminosity L(Lya) = (2.0 + 0.1) x 10 43 erg s _1 . 
The line equivalent width is W Q (Lya) = (11 ± 0.2 ± 2.5) A, 
(switching sign convention, compared to Table|2] to indicate equiv- 
alent widths of emission lines as positive and those of absorption 
lines as negative); here we have quoted separately the random and 
systematic errors respectively to emphasize that the main uncer- 
tainty in the value of W / o(Lya) arises from the placement of the 
underlying stellar continuum. Although the net profile is in emis- 
sion, the value of Wo (Lya) we measure is below the threshold 
Wo (Lya) ^ 20 A generally adopted to define the so-called Lya 
emitters (e.g. Hu et al. 1998; Rhoads et al. 2000; Shapley et al. 
2003). Nevertheless, the profile exhibits the broad characteristics 
that are typical of Lya emission from high-z galaxies: an abrupt 
blue edge and an extended red wing (e.g. Tapken et al. 2007); with 
the high resolution and S/N of our data we can now examine its 
characteristics in greater detail than is normally possible. 

6.2 Comparison to radiative transfer models 

The Lya emission line emerging from star-forming galaxies has 
been the subject of many studies over the last forty years, since 
the seminal paper by Partridge & Peebles (1967). Among more re- 
cent investigations, those by the Geneva group (see, for example, 
Schaerer 2007 which also includes a comprehensive set of refer- 
ences to earlier work) have focussed in particular on the radiation 
transfer of Lya photons in an expanding medium, in an attempt to 
reproduce the variety of Lya profiles seen in high-redshift galaxies, 
where large-scale outflows are the norm. Having established the ex- 
istence of such outflows in the Cosmic Horseshoe (Section [5TTl l, we 
can use the models by Verhamme et al. (2006) to interpret the Lya 
profile in Figure|7] 

The basic idea is that the Lya photons emitted from a cen- 
tral region of star-formation, where the early-type stars and their 
H II regions are located, suffer multiple scatterings before they can 
escape the nebula. Such scatterings not only redistribute the pho- 
tons in frequency, but also convert a fraction of them into infrared 
photons if any dust is present. If the scattering medium is expand- 
ing, the net effect is a reduction of the Lya luminosity (measured 
by the fraction, / csc , of photons which manage to escape) and an 
overall shift of the wavelength distribution of the escaping photons 
to longer wavelengths, as blue-shifted photons are absorbed by the 
Hi gas in front of the stars. The profile of the emergent Lya line 
(for a given star formation rate which determines the intrinsic Lya 
luminosity) depends on a combination of parameters, principally 
the expansion velocity of the medium, i> oxp ; its internal velocity 
dispersion, normally measured by the parameter b = \/2a\ the col- 
umn density of neutral hydrogen, iV(Hl); and the column of dust, 
as measured by the colour excess E(B — V). 

The Lya line in the Cosmic Horseshoe matches well some 
of the characteristics of the theoretical profiles generated by the 
Verhamme et al. (2006) models for the simplest geometrical con- 
figuration: a spherical shell expanding with uniform velocity. In- 
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deed, as can be appreciated from FigureQwhere the Lya emission 
line and the Si II A1526.7070 absorption line are plotted on a com- 
mon velocity scale, the red wing of the former almost mirrors the 
blue wing of the latter, both extending to |«| ~ 700-800 km s _1 . 
A prediction of the models is that there should be a maximum 
in the Lya emission at a positive velocity (relative to the stars) 
v — — 2ii exp . In our case, the velocity of the longer- wavelength 
peak, v — +275 km s -1 is indeed at about twice = 146 kms -1 
which we measured in Section [5"Tl for the gas with the highest opti- 
cal depth in absorption. Referring to Figure 12 of Verhamme et al. 
(2006), these are photons which have undergone just one backscat- 
tering from the receding portion of the expanding shell into our 
line of sight — with this frequency shift they can travel unabsorbed 
through the approaching (i.e. blue-shifted) part of the shell. The 
long tail to v ~ 700 km s -1 is made up of photons which have un- 
dergone two or more backscatterings, while the narrow component 
centred at v = +115 km s -1 consists of photons which have been 
multiply scattered from the approaching part of the shell. 

The balance between the photons escaping from the approach- 
ing and receding portions of the shell — reflected in the relative flux 
of the two emission peaks — is a function of the column density of 
neutral hydrogen through the shell: as iV(H I) increases, the more 
redshifted peak of emission dominates over the lower redshift one. 
Out of the illustrative examples considered by Verhamme et al. 
(2006) in their Figure 16, the case with iV(Hl)= 7 x 10 19 cm -2 
matches our observed profile closely. 

While this correspondence between simulated and observed 
profiles is impressive, other aspects of the models are less satisfac- 
tory. One of the difficulties in relating the Verhamme et al. (2006) 
models to our observations of the Cosmic Horseshoe stems from 
the fact that the absorption profiles of the interstellar lines do not fit 
the simple picture of an expanding shell with a well-defined veloc- 
ity. With absorption extending over a thousand km s~\ and with- 
out any information which would allow us to relate the velocity 
of the gas to its location within the galaxy, there is no obvious 
justification for equating the velocity of the gas with the highest 
optical depth, v — —146 km s -1 , to the expansion velocity of 
the idealised shell. A further difficulty is that in the Verhamme 
et al. (2006) models there is also a strong dependence of the line 
profile on the velocity dispersion of the gas within the expand- 
ing shell. The clear distinction between the two redshifted peaks 
of emission requires not only low column densities of neutral gas 
[iV(H I) ^ 5x 10 20 cm -2 ], but also relatively quiescent gas within 
the expanding shell (|u e x P |/& ~ 5). With the values we measure 
from the interstellar lines in the Horseshoe, \v cxp \/b = 145/105 ~ 
1, the profiles predicted by Verhamme et al. (2006) would look 
drastically different (see their Figure 15) and the two emission com- 
ponents would merge into one. The problem is exacerbated by the 
presence of gas at velocities as high as v ~ -800 km s -1 (see 
Figure|4]l which is not included in the models. A final point is that 
the models considered by Verhamme et al. (2006) generally assume 
unity covering factor of the central source of Lya photons by the 
outflowing gas, whereas in the Cosmic Horseshoe the interstellar 
absorption lines appear to cover only ~ 60% of the stellar UV con- 
tinuum (Section r5.2t . 

Presumably, the real geometrical configuration of stars, Hll 
regions (which recall are at the same velocity as the stars — see Ta- 
bleQJ, and outflowing ISM is more complicated than in the simple 
picture considered (for general illustrative purposes) by Verhamme 
et al. (2006). Treatments of Lya radiation transfer in more com- 
plex situations have been examined in the literature, but generally 
with fewer specific realizations that can be compared directly to 
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Figure 7. Ly« emission and Sill absorption lines plotted on a common 
velocity scale. The former mirrors the latter in velocity space. 

our observations of the Lya line in the Cosmic Horseshoe. For ex- 
ample, Hansen & Oh (2006) considered the effects of a multiphase 
gas outflow, either in a shell with holes or in an ensemble of gas 
clumps. It is interesting that in their simulations the former sce- 
nario can still give a double-peak Lya emission line, whereas the 
latter tends to generate more amorphous profiles. A further point 
of note is that in such multiphase outflows Lya photons can escape 
even when the intervening column density of neutral hydrogen is 
high. It would obviously be of considerable interest to test whether 
a model tailored to the observed properties of the Cosmic Horse- 
shoe can reproduce its Lya profile. 



6.3 Escape fraction of Lya photons 

Before concluding this section, we calculate the escape fraction of 
Lya photons by comparing the observed Lya luminosity with that 
of the UV continuum. In carrying out such a comparison, we con- 
vert the luminosities to star formation rates (SFR) using the conver- 
sions proposed by Kennicutt (1998), divided by a factor of 1.8 to 
account for the lower proportion of low mass stars in the Chabrier 
(2003) IMF relative to the standard Salpeter (1955) IMF adopted 
by Kennicutt. We also correct for gravitational lensing magnifica- 
tion, although our conclusion is not sensitive to this factor, as we 
are comparing luminosities over the same regions of the Cosmic 
Horseshoe. In the lensing model by Dye et al. (2008), the magni- 
fication factor of the whole Einstein ring is 24; in Section [3] we 
calculated that the ESI slit captured a fraction 0.25 ± 0.05 of the 
total light. In the following, we do not carry the uncertainty in this 
estimate through the calculation, as it is the same for the Lya line 
and the UV continuum, and we adopt a magnification factor of 
24 x 0.25 = 6 for the portion of the Cosmic Horseshoe recorded 
with ESI. The comparison between SFR(Lya) and SFR(UV) is sen- 
sitive to the reddening correction, since the two estimates are based 
on features at different wavelengths. We adopt E(B — V) = 0.15, 
as determined for the same two knots of the Horseshoe by Hainline 
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et al. (2009 — see their discussion of the derivation of this value), 
together with the reddening curve of Calzetti et al. (2000). 

From our ESI spectrum we measure a UV continuum lumi- 
nosity at 1700 A L„(1700A) = 1.1 x 10 30 erg s" 1 Hz" 1 which, 
using Kennicutt's (1998) conversion 

SFR(UV)(M yr' 1 ) = 1.4 x 10~ 28 L„ (erg s' 1 Hz" 1 ), (4) 

gives 

SFR(UV) = 159 x 1/1.8 x 3.8 x 1/6 = 56M Q yr" 1 (5) 

with correction factors for the IMF, reddening, and magnification 
respectively^ 

For comparison, Hainline et al. (2009) deduced SFR(Ha) = 
113 ± 17Mq yr _1 from the Ha luminosity measured with NIR- 
SPEC, but the comparison is complicated by the uncertainties in the 
respective flux calibrations. In any case, in galaxies at z = 2 — 3, 
there is an inherent dispersion (of about this magnitude) in star for- 
mation measurements from Ha emission and UV continuum (e.g. 
Pettini et al. 2001; Erb et al. 2006c) which can have a number of 
causes, apart from random errors. First, emission lines and UV con- 
tinuum do not sample the same stellar populations: the luminosity 
of Ha is due primarily to the most massive stars, with shorter life- 
times than the wider range of stellar masses (M £ 5Mq) whose 
integrated light makes up the continuum at 1700 A. Second, the 
conversion from Ha (and Lya) luminosity to SFR implicitly as- 
sumes that all the Lyman continuum (LyC) photons are absorbed 
within the H II region and reprocessed into emission lines — a sit- 
uation which is often referred to as a 'radiation-bounded nebula'. 
Of course, if a fraction of LyC photons escape unabsorbed from the 
H II region, SFR(Ha) and SFR(Lya) will be systematically lower 
than SFR(UV). 

Turning to Lya, we can use Kennicutt's (1998) calibration of 
SFR(Ha) and case B recombination to deduce 

SFR(Lya)(M Q yr _1 ) = 9.1 x 10" 43 L(Lya) (erg s" 1 ) (6) 

which, together with our measured L(Lya) = 2.0 x 10 43 erg s" 1 , 
gives 

SFR(Lya) = 18 x 1/1.8 x 5.3 x 1/6 = 8.9M yr" 1 (7) 

where, as in eq. {5}, the correction factors account for the IMF, dust 
extinction, and magnification respectively. 

By comparing the three estimates of SFR, we deduce ~ 
0.16-0.08, depending on whether we compare SFR(Lya) to 
SFR(UV) or SFR(Ha), respectively. The former value is inde- 
pendent of the absolute flux calibration but assumes that, intrinsi- 
cally (i.e. before the Lya line is quenched by resonant scattering), 
SFR(UV) = SFR(Lya). Conversely, the latter value assumes no er- 
rors in the relative flux calibrations of the NIRSPEC and ESI data, 
and in the calculation of the fractions of the total light from the 
Horseshoe captured by the two spectrographs. However, either esti- 
mate of fc/ c a is plausible, given the wide range found in other Lya 
emitters (e.g. Verhamme et al. 2008; Pentericci et al. 2008). It must 

1 Recently, Leitherer (2008) re-examined Kennicutt's calibrations as given 
here in eqs. (4) and {6) in the light of new Starburst99 models which in- 
clude the effects of stellar rotation. While these models are still at an ex- 
ploratory stage, the indications are that, at solar metallicities, they may re- 
sult in ~ 20% reductions of the values of SFR deduced from both the UV 
continuum and Ha. At sub-solar metallicities, however, SFR(Ha) is fur- 
ther reduced relative to SFR(UV), reflecting the harder ionizing spectrum 
of metal-poor massive stars. Such effects may partly explain why we deduce 
SFR(Hq) > SFR(UV) in the Cosmic Horseshoe. 
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Figure 8. After scaling the flux in Aperture 2 (see Figure [TJ by a factor of 
0.84, the Lya profiles recorded from the two bright knots of emission in the 
Cosmic Horseshoe are indistinguishable from each other at the S/N of the 
present observations. The corresponding error spectra are also plotted (thin 
line near the zero level). 

also be remembered that these values of /^f apply to our view- 
ing angle and to the area of the galaxy covered by the entrance slit 
of the spectrographs; the escape fraction could be higher in other 
directions (Neufeld 1991) and over larger areas (Saito et al. 2006). 

Concluding this section, models of Lya emission resonantly 
scattered by an expanding medium need further improvement in or- 
der to reproduce simultaneously the observed properties of the Lya 
emission line and of the interstellar absorption lines in the Cosmic 
Horseshoe. Nevertheless, interpreting the Lya line in terms of these 
models has given us some insight into the physical conditions of the 
gas. The main inference one would draw from the models of Ver- 
hamme et al. (2006) is that the column densities of neutral hydrogen 
and dust need to be relatively low — iV(H I) — 7 x 10 19 cm -2 and 
E(B — V) ~ 0.1 — in order for the line to exhibit two distinct emis- 
sion peaks at the velocities observed. On the other hand, recall that 
in Section POl we concluded that log iV(H l)/cm~ 2 ~ 20.8 (/ the 
first ions arise primarily in H I gas, as is the case in the Milky Way, 
damped Lya systems, and the gravitationally lensed LBG cB58. 
The order of magnitude difference between these two estimates of 
the column density of gas can be reconciled if most of the gas in 
front of the stars in the Cosmic Horseshoe is ionized, with a neu- 
tral fraction of only about 10%. The profile of the O I A1302.1685 
absorption line, which traces exclusively neutral gas, is consistent 
with this conclusion. 

Given the sensitivity of the profile of the Lya line to the geo- 
metrical and physical properties of the ambient interstellar medium, 
one may expect to see variations between even slightly different 
sightlines through the galaxy. In principle, the lensing magnifi- 
cation of the source offers a good opportunity to recognize such 
small-scale inhomogeneities. And yet no such variations are seen 
in the Cosmic Horseshoe. Recall (Section [3] and Figure [TJ that we 
obtained separate spectra for the two brightest knots in the Einstein 
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Table 4. C in] Emission Lines 



Line 


z 












(kms" 1 ) 


(10~ 18 erg s" 1 cm" 2 ) 


(10 41 ergs- 1 ) 


1906.683 


2.38103 ±0.00014 


62 ± 7 


9.9 ±1.3 


1.29 ±0.16 


1908.734 


2.38127 ±0.00017 


62 ± 7 


8.8 ±1.3 


1.14 ±0.16 



a After correcting for the instrumental resolution (cr instr = 32 km s -1 ), and constrained to 
be the same for both lines. 

b Corrected for reddening and lensing magnification. 



ring, which we labelled 'Aperture 1' and 'Aperture 2'. After scaling 
the latter by a factor of 0.84 (derived empirically and found to be 
the same for the Lya line and the UV continuum), the profiles of 
the Lya line in the two apertures are identical within the limits of 
the noise, as can be appreciated from inspection of Figure [8] Evi- 
dently, the same region of the source is lensed into the two knots 
of the Cosmic Horseshoe. This is all the more puzzling, given that 
Hainline et al. (2009) found that the ratio [N II] A6584/Ha differs 
by almost a factor of 2 (at the ~ 4a significance level) between the 
two Apertures. The [N II]/Hq ratio responds to changes in metallic- 
ity and ionization parameter (e.g. Pettini & Pagel 2004; Kobulnicky 
& Kewley 2004) which presumably are not the same in the regions 
of the galaxy lensed into Apertures 1 and 2 within the Einstein ring 
(knots 'A' and 'D' of Belokurov et al, 2007). Evidently, such varia- 
tions do not have an impact on the Lya emission morphology. 



7 OTHER EMISSION LINES 

The nebular C III] AA1906. 683, 1908.734 doublet is clearly re- 
solved in our ESI spectrum of the Cosmic Horseshoe (see Figure^. 
By fitting Gaussian profiles to the two lines, we deduced the pa- 
rameters listed in Table [4] The mean redshift, zq mj — 2.38115, 
is in excellent agreement with the redshift of the OB stars, as 
expected (Section |4. IK and differs by only —7km s _1 from the 
weighted mean (zhu) = 2.38123 of the two apertures (see Fig- 
ure [TJ measured separately by Hainline et al. (2009) from their 
NIRSPEC spectra. The good match attests to the mutual consis- 
tency of the ESI and NIRSPEC wavelength scales. The line widths 
are rather uncertain due to the limited S/N; with the prior condi- 
tion that the width should be the same for both lines, we obtain 
°"c in] = 62 ± 7km s~ (this is the value used in the fits shown in 
Figure [9}, in very good agreement — again, as expected — with the 
weighted mean (au a ) = 65km s~ x of the two apertures reported 
by Hainline et al. (2009). 

As is well known, the ratio of these two C III] lines is a func- 
tion of the electron density, varying from F( 1906) /F( 1908) = 
1.5 to ~ 0.8 in the range n(e) = 100-30 000 cm -3 . Our measured 
F(1906)/F(1908) = 1.1 ± 0.2 implies very high densities, in the 
range n(e) ~ 5 000-25 000 cm -3 , which however may not be un- 
usual for starbursts at z — 2 — 3 (Brinchmann et al. 2008; Liu et 
al. 2008; Hainline et al. 2009). 

We searched for nebular O III] AA1660.809, 1666.150 emis- 
sion, but these lines are below the detection limit of our data. 
This is not surprising, given that the O III] doublet is weaker than 
C III] by a factor of ~ 7 in the LBG composite spectrum of Shap- 
ley et al. (2003), and the Cm] lines in the Cosmic Horseshoe 
are only detected at the ~ 7a significance level (Table |4). Any 
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Figure 9. Histogram: Portion of the ESI spectrum of the Cosmic Horseshoe 
in the region of the C III] AA1906.683, 1908.734 nebular emission lines, 
after subtracting the stellar continuum. Blue continuous line: Gaussian fits 
to the emission lines; see Tableland Section[7]for details. 

He II A1640.418 emission from Wolf-Rayet stars is also too weak 
to be identified with certainty in our spectrum. 



8 DISCUSSION 

We now consider the results of the previous sections in the light of 
some of the questions raised in the Introduction. 

8.1 Metallicity estimates 

With available instrumentation, only strongly lensed galaxies offer 
the means to cross-check at high redshift the metallicity of early- 
type stars with that of the H II regions that surround them. We nat- 
urally expect the two measures to be the same, since the stars pre- 
sumably formed very recently out of the gas which they now ion- 
ize; thus, this is essentially a consistency check between different 
metallicity indicators, all the more important given the known sys- 
tematic offsets between different emission line measures of chemi- 
cal abundance in H II regions (see, for example, Kewley & Ellison 
2008). Such offsets in turn impact on the interpretation of chemi- 
cal patterns and of trends, such as the mass-metallicity relation and 
its evolution with redshift, which can give us insights into galactic 
chemical evolution. 

In Table [5] we have collected available determinations of el- 
ement abundances in the Cosmic Horseshoe, bringing together the 
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Table 5. Metallicity Comparison 



Method 


Element(s) 


z/z* 


Comments 


R23 


O 


1.5 


H II regions' 5 


N2 


o 


0.5 


H II regions' 


03N2 





0.5 


H II regions' 


1425 


C, Si, Fe 


0.5 


Photospheric, OB stars c 


1978 


Fe 




Photospheric, B stars c 


Civ 


C, N, O, Fe 


~ 0.6 


Stellar wind, O stars d 



a Abundance relative to solar (on a linear scale), using the 
compilation of solar abundances by Asplund et al. (2005). 
b As reported by Hainline et al. (2009). 
c This work (Section [4.2> . 
d This work (Sectionl4~3l 



stellar data from Section|4]with the nebular measurements by Hain- 
line et al. (2009). The first three entries in the Table give the oxy- 
gen abundance deduced from strong emission line indices, respec- 
tively the R23 index of Pagel et al. (1979), for which Hainline et 
al. (2009) used the calibration by Tremonti et al. (2004) based on 
observations of tens of thousands of nearby galaxies in the SDSS, 
and the N2 and 03iV2 indices of Pettini & Pagel (2004). The last 
three entries in the Table are from our analysis in this paper of UV 
spectral features from OB stars in the Horseshoe. 

It is difficult to quote an error for each of the entries in Table[5j 
a realistic estimate should take into account the S/N of the data, 
the scatter in the adopted index calibration and the systematic bias 
which may affect it. However, most estimates evidently converge 
towards a metallicity of about a half solar (which, incidentally, is 
quite typical of star-forming galaxies at these redshifts — see Erb 
et al. 2006a). It is reassuring that stellar and nebular abundances 
are generally in good mutual agreement. The only deviant measure 
appears to be from the R23 method which is known, however, to 
overestimate the oxygen abundance in the near-solar regime, par- 
ticularly when the calibration by Tremonti et al. (2004) is used (e.g. 
Kennicutt, Bresolin & Garnett 2003; Kewley & Ellison 2008). 

Less satisfactory aspects of the comparison in Table[5]are: (i) 
the poorly understood failure of the 1978 index of Rix et al. (2004) 
to provide a meaningful estimate of metallicity in this particular 
case; and (ii) the current limitations in the use of the wind lines 
as accurate abundance diagnostics. Point (i) can only be addressed 
with a larger set of high quality UV spectra of high redshift galax- 
ies. Until it is resolved, however, it may be unwise to use this index 
to investigate the possibility of a differential enrichment of Fe-peak 
elements compared to the products of Type-II supernovae (Pettini 
et al. 2002; Halliday et al. 2008). Concerning point (ii), it is frustrat- 
ing that while the strong wind lines are in principle one of the most 
easily measured abundance diagnostics — at least in data of suffi- 
cient resolution to resolve stellar and interstellar components — we 
still lack a comprehensive library of empirical ultraviolet spectra of 
OB stars of different metallicities to realise their full potential. The 
HST survey of the Magellanic Clouds by Leitherer et al. (2001) 
went some way towards remedying this situation, but their sam- 
pling of the upper H-R diagram is still too sparse to assemble sepa- 
rate sets of Large and Small Magellanic Cloud stars. The resulting 
hybrid library, obtained by combining all the available spectra into 
one set, can only give an approximate measure of metallicity, given 
the factor of ~ 2 difference in the oxygen abundance of the two 
Clouds. 



8.2 Escape of ionizing photons 

It is still unclear what determines /<^c C , the fraction of hydrogen 
ionizing photons which escape from star-forming galaxies into the 
intergalactic medium. Direct detection of these Lyman continuum 
photons has proved problematic until recently (e.g. Shapley et al. 
2006; Iwata et al. 2009), and yet f^ c must have been large at 
very high redshifts for the Universe to be reionized by the star- 
formation activity thought to have taken place at z > 6 (e.g. Bolton 
& Haehnelt 2007; Ryan-Weber et al. 2009). 

Observations of gravitationally lensed galaxies may offer in- 
sights into the factors that control f^ c . By fully resolving the in- 
terstellar absorption lines in our ESI spectrum of the Cosmic Horse- 
shoe, we reached the conclusion that the interstellar gas only cov- 
ers ~ 60% of the stellar UV light, as viewed from Earth. This is a 
promising prerequisite for large values of /<^c C , making the Horse- 
shoe a high priority candidate in searches for LyC emission from 
high-z galaxies. On the other hand, if ~ 40% of the photons from 
the stars and surrounding H II regions really had a clear path out of 
the galaxy, we may have expected to see a strong and narrow Lya 
emission line centred at zhh, whereas no such feature is present in 
our spectrum. Possibly, some 40% of the stars are located behind 
neutral gas of too low a column density to give discernible absorp- 
tion in the metal lines, and yet capable of scattering most of the 
Lya photons out of the line of sight (Hansen & Oh 2006). Whether 
such gas would be optically thick to LyC photons remains to be 
established. 

Another necessary condition for the escape of LyC photons is 
a weaker Ha emission line than expected on the basis of the UV 
continuum luminosity and reddening (admittedly in the idealised 
case of continuous star formation at a constant rate). Such a dispar- 
ity would arise from a 'matter-bounded nebula', where not all of 
the LyC photons emitted by the stars are absorbed and reprocessed 
within the Hll region. As discussed in Section 1631 in the Cosmic 
Horseshoe SFR(Ha) -jt SFR(UV), although a number of different 
factors, apart from leakage of LyC photons, can affect this com- 
parison. In conclusion, it would definitely be worthwhile to search 
for LyC emission in the Cosmic Horseshoe at wavelengths below 
3085 A (the redshifted value of the Lyman limit at z = 2.38115). 
Furthermore, it would be of interest to check for partial covering 
of the stars by the foreground interstellar medium in more galaxies 
among the newly discovered strongly lensed sources, ideally in- 
cluding galaxies with a range of Lya equivalent widths. With such 
data in hand, we should be in a better position to understand the 
conditions that determine the escape fraction of LyC photons and 
the morphology of the Lya line. 



8.3 Comparison with MS 1512-cB58 

One of the motivations for the present study was to establish how 
typical are the properties of the galaxy MS 1512-cB58, the only 
previous case where the gravitational lensing boost was sufficient to 
allow a detailed look at the spectrum of a high-redshift star-forming 
galaxy. While we have now doubled the 'sample' of high-z galax- 
ies with good-quality ESI spectra, it would clearly be premature 
to draw general conclusions on the basis of just two objects. Nev- 
ertheless, one cannot help but being struck by how closely these 
two galaxies resemble each other in many of their properties. They 
are very similar in their overall metallicity and probably in their 
detailed chemical composition, indicating that they have reached 
comparable stages in the conversion of their gas reservoirs into 
stars. Their young stellar populations are largely indistinguishable 
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in their UV spectral features, especially the P-Cygni lines which 
are sensitive indicators of the mode of massive star formation and 
of differential reddening among stars at the upper end of the IMF. 
Most remarkable perhaps is the unexpected similarity in the proper- 
ties of their outflowing interstellar media, with velocities spanning 
in both cases ~ 1000 km s _1 , from —800 to +250 km s _1 relative 
to the stars, and involving comparable column densities of gas. 

Presumably, it is just a coincidence that the first two galax- 
ies at z = 2-3 whose UV spectra have been put under the 'mi- 
croscope' of gravitational lensing have also turned out to be sim- 
ilar in so many respects. After all, we now know that the Cosmic 
Horseshoe and MS 15 12-cB58 are also forming stars at comparable 
rates of ~ 50-100 Mq yr _1 , and have similar dynamical masses 
of Afvir ~ 1 x 10 10 M S (Hainline et al. 2009 for these measure- 
ments in the Cosmic Horseshoe and Teplitz et al. 2000 for cB58). 
Future studies of other gravitationally lensed objects will arguably 
show a wider range of properties of star-forming galaxies at these 
redshifts. 

What is interesting then, given the underlying similarity of the 
Horseshoe and cB58, is to consider the aspects in which they dif- 
fer, the incomplete coverage of the stellar light by the interstellar 
absorption lines and the strikingly different Lya morphologies. It is 
these two factors alone that would lead an observer to 'classify' the 
two galaxies differently. For example, in the rough breakdown of 
LBGs by Shapley et al. (2003) according to their UV spectral prop- 
erties, cB58 would fall in the quartile with the highest interstellar 
absorption line equivalent widths and reddening, while the Cosmic 
Horseshoe would be in the quartile with the second highest Lya 
emission equivalent width. And yet, from what we have uncovered 
in the present study, it may be the case that these differences are 
perhaps simply due to orientation effects — viewing two intrinsi- 
cally similar galaxies along differing sightlines through their inter- 
stellar media, one (the Horseshoe) more rarefied and highly ionized 
than the other (cB58). With future studies of other strongly lensed 
galaxies we shall hopefully learn how to distinguish intrinsic dif- 
ferences from those caused by second-order effects, of which the 
direction from which we view the central starburst may be one ex- 
ample. For the moment, however, we must remain cautious of over- 
interpreting differences between different 'classes' of high-redshift 
galaxies, especially if such distinctions are based on the morphol- 
ogy and strength of the Lya line, whose appearance is so sensitive 
to a variety of different parameters (e.g. Pentericci et al. 2009). 



9 SUMMARY AND CONCLUSIONS 

We have presented ESI observations of the rest-frame ultraviolet 
spectrum of a star-forming galaxy at z — 2.38115, lensed by a mas- 
sive foreground galaxy into a nearly complete Einstein ring dubbed 
the 'Cosmic Horseshoe'. The high gravitational magnification af- 
fords the rare opportunity of recording the spectrum of this member 
of the z — 2-3 galaxy population at higher resolution and S/N than 
otherwise achievable, and thereby provides us with a close-up look 
of a galaxy during the epoch when cosmic star-formation activity 
was at its peak. Our main findings are as follows. 

(i) The Cosmic Horseshoe shares many of its properties with 
the population of 'BX' galaxies selected with the colour criteria of 
Steidel et al. (2004). Its metallicity Z ~ 0.5Zq, dynamical mass 
M vir ~ lxlO lo M , and reddening E(B- V) = 0.15 are all typ- 
ical of those galaxies (Erb et al. 2006a,b,c). With a star-formation 
rate SFR~ lOOAf© yr _1 , the Horseshoe is among the most lumi- 
nous galaxies in the BX sample. 



(ii) Generally, there is good agreement between different 
metallicity indicators based, respectively, on ratios of strong neb- 
ular emission lines, on blends of stellar photospheric lines, and on 
P-Cygni lines from the most luminous OB stars. The R23 index of 
Pagel et al. (1979) seems to overpredict the oxygen abundance by a 
factor of ~ 3, when calibrated with the prescription by Tremonti 
et al. (2004). Of the stellar measures, the 1978 index of Rix et 
al. (2004) does not fit the data well due to excess absorption not 
present in the models and whose identity remains unclear. 

(iii) A continuous mode of star formation with a Salpeter slope 
for stars with masses 100 ^ M ^ 5 Mq gives a good representa- 
tion of the UV spectrum; we see no evidence for a departure from 
a Salpeter IMF, nor differential reddening among massive stars, in 
the contrast of the P-Cygni lines over the underlying UV contin- 
uum. 

(iv) The interstellar absorption lines are broad, extending over 
a velocity interval of ~ 1000 km s~\ from -800 to +250 km s~ x 
relative to the redshift of the stars and H II regions. The gas with 
the highest optical depths is outflowing from the regions of star 
formation with a speed of ~ 150 km s , and only covers ~ 60% 
of the stellar continuum light. 

(v) The Lya line is in emission and shares many of the prop- 
erties of the so-called Lya emitters, although its equivalent width 
is a factor of ~ 2 smaller than the lower limit Wo (Lya) ^ 20 A 
commonly adopted to define this subset of galaxies. The resolved 
line profile matches well those computed with models of reso- 
nantly scattered Lya photons in an expanding medium, although 
such models find it difficult to account for the fact that the velocity 
dispersion of the gas is comparable to its bulk expansion velocity, 
and for the existence of gas extending over hundreds of km s^ 1 in 
velocity space. Some ~ 10-15% of the Lya photons produced es- 
cape the galaxy, the remainder presumably being absorbed by dust 
and recycled into infrared photons, or scattered over a larger area 
than that covered by our narrow-slit observations. It remains to be 
established what fraction of the Lyman continuum photons emitted 
by the stars leak into the IGM. 

(vi) Overall, many of the physical properties of the stars and 
interstellar medium of the Cosmic Horseshoe are similar to those 
of the only other galaxy previously studied at comparable resolu- 
tion and S/N, MS 1512-cB58 at z = 2.7276. This may not be sur- 
prising, given that both galaxies are fairly typical examples of the 
star-forming galaxy population at these redshifts. The fundamental 
similarities between the two objects highlight the fact that the as- 
pects in which they differ, particularly the very different morpholo- 
gies of their Lya lines, may be due to superficial reasons, such as 
differing viewing angles, rather than to more deeply rooted causes. 

In the near future it should be possible to make further 
progress on all of these topics with the increasing numbers of 
strongly lensed high-redshift galaxies being discovered by the 
SDSS and other lensing surveys. 
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